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mltlvarlabla  faedback  coocrol  arttaa  la 
dengnad  for  a  subaaralbla.  Tha  control  varlablaa 
ara  tha  bow,  ruddar,  and  dlffaraneial  team  control 
aurfacaa;  thaaa  ara  dynaalcally  coordlnatad  to  aa 
to  cauaa  tha  vahlcla  to  follow  Indapandant  and 
ataultanaoua  eowwandad  changaa  In  yaw  rata,  dapth 
rata,  pitch  attltuda,  and  roll  angla.  Two  daalgna 
wara  aaaluatad  utlng  a  nonllnaar  aubwartlbla  alau- 
latlon.  Ona  utad  all  four  control  varlablaa  to  that 
actlva  roll  control  waa  poaalbla.  Tha  othar  uaad 
only  throa  control  varlablaa,  and  actlva  roll 
control  waa  not  taployad.  Both  faadback  ayatawa 
wara  daalgnad  ualng  the  Llnaar  Quadratic  Gauaalan 
with  Loop  Tranafar  Racovary  ^(LTRT^  daalgn 
aathodology  to  at  to  naat  alallar  daalgn  apaclflca- 
tlona  In  tha  fraquaney  doaaln.  Both  tha  llnaarlaad 
aodala,  and  tha  non-llnaar  alaailatlon  hava  ahown 
that  actlva  roll  control  ylaldt  a  vary  algnlfleant 
laprovaaant  In  aubwartibla  parfonanea.  Actlva 
roll  control  alnlalaad  unwanted  dapth  changaa  In 
difficult  cowaandad  trajectory  acanarloa. 

^  I  ,  / 

INTRODU^ION 

Thla  paper  auaoarlzea  the  raaulta  ot  a  feaalbt- 
llty  atudy,  (1),  related  to  the  daalgn  of  multi¬ 
variable  feedback  control  ayateaa  for  deeply 
aubmarged  aubnaralblea.  Preaant  aubmeralblaa  have 
minimal  automatic  control  capabllltlea.  Reaearch 
and  development  af forte  ara  underway  to  examlna  the 
banaflta  of  eloaad  loop  control  for  maneuvering 
aubnaralblaa.  A  key  laaua  that  haa  received  little 
attention  la:  what  are  the  banaflta  of  providing 
for  Independent  actlva  roll  control  of  a  aubmer- 
albla  by  modification  of  tha  atarn  control 
aurfacaa?  The  objectlva  of  our  feaalblllty  atudy 
waa  to  examine  thla  Important  laaua. 


1.  Reaearch  conducted  at  the  NIT  Laboratory  for 
Information  and  Declalon  Syatema,  wl(h  aupport  pro¬ 
vided  by  the  Office  of  Naval  taaaarch  under 
Contract  ONR/N00014-82-K-OS82(NR60b-003). 

2.  Naval  Sea  Syatema  CoiuMnd,  Code  3SU3t, 
Waahlngton,  D.  C.  20362. 

3.  Department  of  Aaronautlca  and  Aatronautlca, 
NIT,  Cambridge,  NA  02139. 

4.  Department  of  Electrical  Engineering  and 
Computer  Science,  NIT,  Cambridge,  NA  02139. 


Tha  medal  of  tha  aubmerelbla  uaad  In  thla  atudy 
uaaa  a  ruddar  and  atam  plana  configuration  of  the 
ao-callad  cruciform  type.  Submaralblaa  uaa  a  cru¬ 
ciform  atera  with  machanlcally  alavad  upper  and 
lower  ruddare,  and  machanlcally  alavad  port  and 
atarboard  atam  planea.  Tha  advantage  of  thla  cm- 
clfom  atam  plana  configuration  la  that  It  allowa 
for  Intuitive  actlona  on  tha  part  of  a  human  opera¬ 
tor  for  daalrad  ahlp  motion.  For  example.  If  It  la 
daalrad  to  rlae  or  diva,  all  the  human  operator  haa 
to  do  la  command  rlae  or  dive  on  tha  atam  planea. 
A  almllar  altuatlon  occura  If  tha  operator  daalrae 
to  turn. 

A  major  drawback  to  thla  atam  configuration  la 
that  chare  la  no  opportunity  to  actively  control 
tha  aubmaralbla  roll  angla.  A  aubmeralbla  haa  a 
natural  candancy  to  roll  by  a  algnlfleant  amount 
during  a  turn.  The  ao-callad  *anap  roll*  of  a 
turning  aubmarlna  la  a  nonllnaar  function  of  tha 
apaad  and  tha  rudder  dlaplaeamant  angla,  and  la 
Induced  by  complex  hydrodynamic  forcaa  and  momenta, 
Tha  detrimental  affaeta  of  non-taro  roll  anglaa  ara 
due  to  the  fact  that  they  reault  In  undaalrabla 
depth  changaa,  and  It  bacomaa  even  more  difficult 
to  turn  while,  at  Che  aaaa  tine,  maintaining  a 
desired  pitch  and/or  depth  rata  command.  Hence,  If 
methoda  are  utilised  Co  reduce  Che  anap  roll,  the 
ability  to  .eaincaln  ordered  depth  during  a  maneuver 
la  graatly  enhanced.  It  la  important  to  keep  In 
mind  Chat  the  large  length  of  a  aubmaralble  (over 
three  hundred  feet)  cranalates  Co  a  very  limited 
tolerance  for  depth  errora  chat  compromlae  tha 
safety  of  the  submersible  and/or  lea  missions. 

By  unslavlng  the  port  and  starboard  planes,  and 
by  commanding  their  daflactlona  Independently,  It 
is  poaalbla  to  actively  control  tha  aubmeralbla 
roll  angle  while  maintaining  pitch  control.  Also, 
the  daf lections  of  the  bow  (falrwater)  plane  can  he 
coordinated  with  Chose  of  tha  stern  planes  to  pro¬ 
vide  additional  control  In  the  longitudinal  axis. 
It  is  obviously  hard  for  a  human  operator  to  com¬ 
mand  coordlnatad  changes  In  four  control  surfaces 
slmultanaoualy .  This  provides  the  motivation  for 
studying  the  automatic  multivariable  feedback 
control  problem  for  maneuvering. 

The  majority  of  previous  controller  designs  for 
submerslblas  used  classical  single-input  single- 
output  (SISO)  techniques,  decoupling  the  longitudi¬ 
nal  and  lateral  axes,  and  using  only  stem  and 
rudder  control.  There  have  been  a  limited  number  of 
multi-input  multi-output  (NINO)  designs  for  full 
scale  submaralblaa  that  can  be  found  In  the  open 
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llt«raCttr«t  la  Ch«M  atudlaa  acclva 
eoncrel  of  cha  roll  angla,  ualng  dlffaranclal  acara 
plana  daflacctona,  waa  not  conatdarad. 

In  Chta  papar  wa  uaa  tha  MIMO  LQG/LTR  daalgti 
nathodology,  (2),  (.3),  Co  daalgn  two  aultlvarlabla 
faadback  conerol  ajracaaa.  Oaa  daalgn  utlllaaa 
ladapandane  eoncrol  of  Cha  pore  and  aCarboard  acam 
planaa  co  achtava  actlaa  roll  control,  whila  the 
other  daalgn  doaa  not.  Tha  two  faadback  apataaa 
are  auppoaad  to  cauaa  cha  aubwaralbla  to  follow 
deaandlng  aanauvarlng  trajaccorlaa  In  raaponaa  to 
alaulcanaoua.  coawaada  In  daalrad  pitch,  depth  rata, 
and  yaw  rata.  The  daalrad  roll  angle  la  aero. 
Slallar  coawand-followlng  apaclflcatlona  ware  poaad 
for  both  daalgna  In  the  frequency  domain.  Tha 
bandwidth  waa  aelacted  ao  chat  axceaalva  aaturaclon 
of  cha  control  aurfacaa  waa  avoided.  In  both  caaaa 
cha  daalgn  aodal  for  cha  aubwaralbla  waa  obtained 
by  llnaarlalng  the  nonlinear  differencial  equaclona 
of  notion  at  high  apaad.  Tha  raaulclng  LQC/LTR 
daalgna  warn  evaluated  ualng  a  full-acale  nonlinear 
alwulatlon  of  tha  aubwaralbla  dynawlca  at  the 
Char lea  Stark  Draper  Laboratory. 

Tha  alwulatlon  reaulta  denonatrate  that  cha  uaa 
of  active  roll  eoncrol  can  have  a  vary  algnlflcanc 
tnpacc  on  the  perforwanca  of  the  aubmeralbla,  even 
In  difficult  Cuming  naneuvera  for  a  huian  opera¬ 
tor.  Active  roll  control  can  reduce,  by  a  leaat  a 
factor  of  two,  undealrabla  depth  excuralona  when 
conparad  to  cha  daalgn  chat  doaa  not  actively 
control  the  anap  roll.  Aiao  va  have  found  that  the 
LQG/LTR  cowpanaacor  la  quite  robuac  co  the  changea 
In  tha  aubwaralbla  dynawlca  during  a  nanauvar, 
where  actual  apeada  can  drop  by  aa  auch  aa  35  per¬ 
cent.  Thua  we  conjecture  that  a  full-envelope 
design  can  be  confidently  carried  out  by  slaply 
gain-scheduling  Che  paraaeCers  of  the  LQC/LTR  com¬ 
pensator  as  a  function  of  aeasured  submersible  for¬ 
ward  speed,  along  the  lines  suggested  by  Lively, 
(4). 

The  remainder  of  this  paper  Is  organized  as 
follows.  Section  2  overviews  the  modeling  Issues 
and  the  definition  of  the  linearized  dynaalcs  used 
to  design  the  feedback  controller.  In  Section  3  we 
present  the  scep-by-scep  procedure  that  we  have 
used  CO  derive  the  multivariable  LQG/LTR 
coapenaacors;  we  augment  the  natural  submersible 
dynawlca  with  Integrators  so  as  co  enforce  zero 
steady-scate  errors  to  constant  comaands  and/or 
disturbances.  The  quality  of  Che  linearized 
designs  la  Judged  by  the  shapes  of  the  singular 
values  of  the  loop  and  closed-loop  transfer  func¬ 
tion  natrlces  In  the  frequency  domain.  We  show  how 
to  select  the  weighting  matrices  associated  with 
Che  LQC/LTR  design  ao  Chat  the  desired  target 
designs,  defined  by  the  singular  values  of  cne 
Kalwan  Filter  loop,  have  Identical  behavior  at  both 
low  and  high  frequencies.  Section  4  presents  and 
Interprets  the  results  of  Che  nonlinear  simulations 
and  contains  a  discussion  of  the  quantlcatlve  ben^ 
fits  of  using  active  roll  control. 

In  closing,  we  stress  chat  the  results  pre¬ 
sented  In  this  paper  represent  only  a  limited 
feasibility  study  on  the  potential  benefits  of 
ualng  multivariable  control  In  general,  and  of 


active  roll  eoncrol  In  particular,  for  full  senla 
subwnrslblas.  The  resultant  LQG/LTR  designs  must  be 
further  refined  and  analyzed  before  they  are  Iwple- 
wanced  In  an  actual  subaaralbla  control  ayatew. 

MODELING 

The  controller  design  procedure  begins  with  Che 
expression  of  Che  equations  of  notion  In  linear 
cine  Invariant  state  apace  fom.  The  nonlinear, 
aultlvarlabla  systaa  chat  represents  the  subaer- 
slble  la  deserllMd  by: 

^x(t)  -  f{x(t),u(t)) 
dt  ■  (1) 

jr,(t)  -  gWt)) 

where: 

£(c}  la  Che  state  vector 
^(t)  la  Che  control  vector 
£(c)  la  Che  output  vector 

The  state  vector  for  the  subaerslble  must 
Include  Che  six  degrees  of  freedow  from  Che  vehicle 
coordinate  ayatew,  Che  desired  position  variables 
CO  locate  Che  subwerslble  with  respect  to  a  fixed 
coordinate  system,  and  the  three  Euler  angles  which 
describe  the  relationship  of  the  motion  of  Che  sub- 
warslble  with  respect  to  the  two  coordinate 
systems.  Figure  1  shows  Che  subaerslble  reference 
systaa  used.  Raferencaa  (9)  co  (12)  provide 
guidance  on  the  developwenc  of  the  subaartne  dyna- 
alca. 

TABLE  I  DEFINITIONS  OF  SUBMERSIBLE  STATES  AND 
CONTROL  VARIABLES 

Submersible  States 


U  •  Xl(t) 

forward  velocity 

(f c/sec) 

V  •  X2(t) 

lateral  velocity 

(f c/sec) 

w  -  X3(t) 

vertical  velocity 

(ft/sec) 

P  •  X4(t) 

roll  rate 

(rad/sec) 

q  •  X5(c) 

pitch  rate 

(rad/sec) 

r  •  xs(c) 

yaw  rate 

(rad/aec) 

♦  •  X7(t) 

roll  angle 

(radians) 

9  -  X8(t) 

pitch  angle 

(radians) 

Control  Variables 


6b  bow/f alrvater  planes  (radians) 

6r  rudder  deflection  (>  left)  (radians) 

6si  port  stern  plane  defection  (radians) 

(■♦•  down) 

6s2  starboard  stern  plane  deflection  (radians) 

(♦  do%in) 
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will  b«  lnv«aclgac«d<  Tha  four  output  aarlablaa 
aalaetad  ara  daptb  rata  i,  paw  rata  i  roll  aogla  d, 
and  pitch  aogla  9« 


With  tha  output  varlablaa  dotamlnad,  and  tha  A 
and  B  natrleaa  ealeulatad  frow  tha  llnaarlaatlonT 
tha  Ttata  apaea  daaerlptlon  of  tha  «odal  la  now 
co^>lata  and  takaa  tha  font 


x(t)  ■  A*(t)  ♦^^(t) 
jr(t)  -  C*(t), 

whara  tha  output  aaetor  ^.it)  la  glvan  bp 

i(t)  -  [dU)  e(t)  Ut)  i(t)f  . 


(2) 


(3) 


rigura  1  Skatch  ahowlng  poaitlaa  dlractlona  of 
axoti  anglaa,  volocltlaa,  forcaa, 
and  woMnta 


Thaaa  nonllnaar  dlffarantlal  aquatlona  can  Chen  ba 
linaarlcad  through  tha  following  fairly  acralght- 
forward  nuaarlcal  cachnlqua*  A  noalnal  point  la 
choaan  for  tha  daalgn  bp  Incagraclng  tha  nonlinear 
aquatlona  of  aoclon  using  a  spaclflad  sec  of  Ini¬ 
tial  conditions*  An  aqulllbrlua  point  Is  found 
chat  corrasponds  to  alnloua  accelerations  for  all 
tha  state  variables  determined  froa  the  Integration 
of  Che  equations  of  aoClon*  The  values  of  Che 
State  variables  at  the  equilibrium  point  then  spe¬ 
cif  p  a  noalnal  point*  about  which  higher  order 
caraa  asp  ba  neglacted.  Proa  chase  results*  a  sat 
of  linear  differential  aquatlona  aap  Chen  be  pr^ 
ducad*  the  A  and  ^  aatrlcas  ealeulatad*  and  a  state 
space  das?rlpclon  of  the  subaerslbla  nodal 
produced. 

The  noalnal  point  chosen  for  this  research 
corresponds  to  a  level  subnerslble  trajectory  at 
high  speed.  The  rudder  deflection,  6r,  cen  be  sec 
at  arbitrary  angles  to  cause  the  submersible  to 
turn  at  different  rates,  end  to  roll  at  different 
angles.  This  accenpcs  to  deteralne  the  open  loop 
sensitivity  of  the  subnerslble  to  roll,  which  has  a 
significant  effect  on  the  depth  of  a  submersible  In 
a  turning  maneuver. 

For  each  noalnal  point  deCeralned,  the 
resulting  linear  nodal  was  validated  by  perturbing 
the  noalnal  point  of  the  linear  and  nonlinear 
nodels*  and  coaparlng  the  tine  histories  of  the 
state  paraaetars.  Provided  tha  perturbations  are 
not  excessive,  the  nonllnaar  nodal  will  return  to 
tha  equilibrium  point.  The  linear  aodel,  however* 
will  not  return  to  Its  equilibrium  point  resulting 
from  the  non-zero  forces  laposed  by  the  control 
forces,  and  the  absence  of  nonllnaar  hydrodynamic 
effects. 


As  nencloned  previously,  there  ara  four 
possible  control  variables.  These  ara  6s j*  6s,* 
6r*  and  6b.  Figure  2  Illustrates  tha  control  suN 
face  configurations  used  In  this  research. 

The  use  of  nodal  analysis  allowed  the  for- 
aulatlon  of  the  prerequisites  necessary  to  pursue 
the  LQC/LTR  design  nechodology.  These  prerequisi¬ 
tes  are  chat  the  open  loop  linear  model  ba  detec¬ 
table  and  stablllzable,  and  chat  the  location  of 
nonalnlaua  phase  zeroa  be  knomi.  For  a  detailed 
description  of  the  analysis  of  the  linear  model* 
refer  to  reference  (1).  The  numerical  values  for 
the  A  and  _B  aatrlcas  are  Included  In  Appendix  A. 
The  7pen-loop  aodel  of  the  linearized  submersible 
dynamics  Is  stable  and  nlnlnum-phasa. 

Tha  final  Issue  associated  with  the  linear 
aodel  Is  chat  of  sealing.  It  Is  now  widely  appre¬ 
ciated  chat  In  multivariable  daslgna  tha  control 
and  output  variables  nusC  be  scaled  so  that 
controls  and  tracking  errors,  expressed  In  dif¬ 
ferent  physlcel  units,  can  be  compared  fairly. 

The  scaled  control  verlables  were  established 
as  follows.  Ue  used  physical  saturation  Halts  on 
Che  control  surfaces:  20  degrees  for  the  bow 
planes,  30  degrees  for  the  rudder,  end  25  degrees 
for  the  stern  planes.  After  transforming  froa 
radians  to  degrees,  the  scaled  controls  were  sec  as 
follows; 

one  unit  of  bow  •  0.667  units  of  rudder 
•  0.8  units  of  stern. 

The  scaled  outputs  were  first  derived  by 
chenglng  radians  to  degrees.  Then,  In  terms  of  the 
performance  of  the  control  systea.  It  was  Judged 
that  a  I  deg/sec  error  In  either  yaw  rate  or  depth 
rate  was  10  tines  more  serious  Chen  a  1  degree 
error  In  roll  or  pitch  engle.  Thus  Che  scaled  ouCv 
puts  were  sec  as  follows: 


The  conperlsons  of  selected  nonlinear  and 
linear  models  and  state  variables  showed  excellent 
correlation,  thus  valldeclng  the  linearized  nodels. 
For  further  details,  refer  to  reference  (1). 

Since  It  Is  desired  to  control  the  submersible 
during  maneuvering  slcuaclona,  a  race  controller 


one  unit  of  pitch  •  one  unit  of  roll 

•  0.1  unit  of  yaw  rate 

•  O.l  unit  of  depth  rate. 

The  resultant  scaled  aodel  wes  used  In  Che 
designs  chat  will  be  described  In  the  sequel. 
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Figure  2  View  froa  Stern  Showing  Rudder 
and  Differential  Stemplanea 


MULTIVARIABLE  CONTROL  STSTEM  DESIGN 
Introduction 

Aa  atated  earlier,  a  controller  la  dealgnad 
ualng  the  LQG/LTR  deaign  aethodologjr*  The  singular 
value  loop  ahaplng  approach  la  used  to  obtain 
daalrable  singular  values  of  the  systea  transfer 
function  aatrlx  to  aeet  the  specifications  of  per- 
foraance  and  robustness  to  plant  uncertainties  and 
aodellng  errors. 

Oontroller  Specifications 

Perforaance  specifications  outlined  In  this 
research  are  not  all  encoapasslng  and  reflect  a 
alapllfled  subaerslble  control  systea.  The  perfo^ 
aance  requlreaents  are  driven  by  the  Intuitive 
engineering  approach  to  obtain  good  coaaand 
following,  good  systea  response,  robustness,  and 
disturbance  rejection.  These  perforaance  requlre¬ 
aents  are  aet  through  loop  ahaplng  techniques. 

Two  performance  requlreaents  are  Imposed  on  the 
controller  design.  First,  the  steady  state  error 
to  step  coaaanda  and  step  disturbances  is  to  be 
zero.  Second,  the  maxlaua  crossover  frequency  Is 
Halted  by  the  ability  of  the  submersible  to 
respond  and  by  the  rate  at  which  the  compensator 
deflects  the  control  surfaces. 


reference  coaaanda  and  the  assaurad  outputs,  and 
was  selected  to  be  In  the  neighborhood  of  0.5 
rad/see.  based  on  desired  actuator  dynaales  and 
saturation  Halts  (1).  Although  not  explicitly 
stated  as  a  perforaance  specification,  froa  the 
perforaance  aspect.  It  la  desirable  to  have  all 
singular  values  cross  over  at  about  the  saaa  fre¬ 
quency.  Also,  on  the  high  frequency  aide,  the 
controller  auat  be  capable  of  rejecting  nolae  and 
be  robust  to  high  frequency  aodellng  errors.  Noise 
sources  generally  originate  froa  the  envlronaent, 
or  froa  the  sensor  Itself.  Sensor  noise  typically 
occurs  at  a  hlghar  frequency  than  the  systea  band¬ 
width  and  should  not  affect  the  dynaales  of  the 
subaerslble  since  ship  eigenvalues  will  typically 
He  In  the  lower  frequency  band. 

The  LQC/LTR  Design  Msthodology 

The  aultlvarlable  LQC/LTR  design  aethodology 
consists  of  four  major  steps  (2),  (3). 

The  first  step  la  the  development  of  a  low  fre¬ 
quency  model  of  the  nominal  plant  and  dateralnatlon 
of  aodellng  uncertainties. 

The  aodellng  uncertainty  In  the  noadnal  node! 
due  to  sensor  noise,  unaodeled  submersible  dyna¬ 
mics,  and  actuator  dynaales.  Is  assiaed  to  be  con¬ 
centrated  at  high  frequencies  (In  this  study, 
higher  chan  10  rad/sec).  Fixing  the  crossover  fre¬ 
quency  of  Che  singular  values  of  the  loop  transfer 
function  matrix  will  detendne  the  significance  of 
the  unaodeled  dynasdes,  and  the  ability  of  Che 
plant  to  aaet  command  following  specif Icaclons. 

The  second  step  of  Che  design  process 
establishes  Che  low  frequency  perforaance  requlrw* 
aencs.  The  state  space  block  dlagraa  of  the  cor 
pensaced  plant  Is  shown  In  Figure  3. 


The  zero  steady  state  requlreaenc  Is  met  by 
placing  Integrators  In  each  of  the  control  char 
nela.  This  Is  necessitated  because  we  plan  to  use 
Che  version  of  Che  LQC/LTR  design  aethodology  which 
breaks  the  loop  at  the  plant  output,  or  equiva¬ 
lently,  at  the  error  channel.  In  this  manner,  the 
Integrators  will  Chen  becone  part  of  the  compen¬ 
sator  which  precedes  the  plant  In  the  feedback 
loop.  Note  chat  the  use  of  Integral  control  In  the 
Input  channels  does  not  prevent  Che  specification 
for  maxlaua  crossover  frequency  froa  being  met. 

The  aaxlaua  crossover  frequency  of  the  compen¬ 
sator  decaralnas  the  rapidity  of  the  control  sur¬ 
face  deflections  based  on  the  error  signals  which 
are  generated  by  the  difference  between  Che 


Figure  3  Block  Dlagraa  of  a  MIMO  Corapensacad  Plant 


r,(s)  •  reference  signal  or  comaand  Input  vector 

_e(s)  ■  error  signal  vector 

u^(s)  -  control  vector  to  the  plant 

^(s)  -  output  vector  of  Che  plant 

^(s)  •  disturbance  vector  at  the  plant  output 

7(s)  >  coapensator  transfer  function  aatrlx 

£(s)  -  augmented  plant  transfer  matrix 

The  transfer  matrix  £(s)  contains  the  noalnal 
low  frequency  model  ^(s)  and  any  augmenting 
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dyaasle*  Gf(»),  tad  1«  daflnad  tha  noalnal  daalgn 
■odalt  Thus 

.£(•)  -  Cp(s)  C,(a)  (4) 

Augasnclng  cha  dynaalcs  of  Cha  aufaMralbla 
concrol  syscaa  noraally  sarvas  a  dual  purposa.  Ona 
la  CO  partially  aodal  tha  actuator  dynaalcs  to  aaka 
tha  aodal  as  aeeuraca  as  posalbla  and  to  achlava 
daslrabla  rolloff  bayond  croasovar  for  atablllCy- 
robuacnass.  Tha  ochar  la  to  Includa  latagrators  to 
causa  tha  coapansator  to  paralc  Cha  subaarslbla  to 
achlava  saro  scaady  atata  arror  to  atap  Inputs  and 
dlsturbancss  (l.a.t  parfacc  coaaand  following  at 
DC).  Tha  'actuator  dynaalca  ara  abova  cha  aazlaua 
axpactad  crossovar  fraquancy,  and  thus  ara 
naglactad.  This  la  parfactly  valid  as  long  as  cha 
rolloff  abova  crossover  la  fast  enough  and  sacla- 
fles  the  robuacnaaa  criteria. 

A  block  dlagraa  of  the  augtMnced  aodel  appears 
In  Figure  4.  It  la  seen  Chat  Che  Incegracora  are 
placed  In  tha  concrol  channels.  The  aathewaclca  of 
Che  augaenced  acacea  are  aanlpulaced  In  such  a  way 
as  Co  provide  a  naans  to  achieve  the  desired  Kalman 
Filter  open  loop  shapes 

i 


Figure  4  Integrators  Placed  In  the  Concrol 
Channel  of  Plant 


crossover,  cha  response  of  Cha  outputs  with  raspaet 
CO  sensor  nolsa  can  be  nlnlnlsad  and  stability- 
robustness  enhanced  If  Cha  aaxlnua  singular  valua 
of  £(s)K(s)  la  snail  with  raspaet  to  unity. 

Combining  cha  above  conditions,  ws  assentlally 
Inposa  high  and  low  fraquancy  barriers  on  cha 
singular  valua  plots  of  the  loop  transfer  function 
matrix  £(s>X(s),  as  shown  In  Flgura  S. 


Figure  5  Plot  of  Desired  Singular  Valua  Shapes 


Tha  high  frequency  barrier  Imposes  a  robuacnaaa 
constraint  on  the  compensator  and  tha  low  frequency 
barrier  Imposes  Che  command  following  and  distur¬ 
bance  rejection  requirements. 


We  define  the  augmentation  dynaalcs  by  ^(s). 
whose  state  space  description  Is 

ap(8)  -  Uc(s)  ;  ^(s)  -  Us 

where  each  oatrlx  la  *^4  x  4 1.  The  augmenting 
Integrator  dynamics  are  added  to  Che  8^"  order 
system  producing  a  12***  order  system.  Note  that 
Che  physical  Input  to  Che  plant  Is  labelled  j^(s) 
CO  distinguish  It  from  the  output  of  the  compen¬ 
sator  U|e(a).  Although  Che  augmentation  dynamics 
^(s)  will  eventually  be  lumped  with  the  compen¬ 
sator,  they  are  kept  separate  until  the  LQC/LTR 
procedure  Is  complete. 


The  third  step  of  Che  design  process  Is  deter¬ 
mining  the  coiq>ensaCor  transfer  function  matrix, 
l((s),  chat  will  provide  the  desired  singular  values 
of  £(s)K(s).  This  step  of  the  process  Is 
appropriately  termed  ’loop  shaping*. 

The  Kalman  Filter  (KF)  methodology  Is  first 
applied  to  Che  nominal  design  model.  This  produces 
a  transfer  matrix  0|(f(s)  that  has  the  desired 
singular  value  loop  shapes;  thus  Is  called  the 
'target  design*.  A  distinction  Is  noted  In  this 
procedure,  however,  because  the  RF  theory  Is 
applied  In  a  specific  manner  which  la  not  to  be 
confused  with  optimal  state  estimation. 


To  determine  the  requirements  of  J((s),  the 
overall  loop  transfer  function  of  the  closed  loop 
system  Is  analysed,  where 

X(s)-l£vG(s)JC(s)  r^d(e)+f_l4C(s)K(s)  f^GlslKCslrls) . 

-  -  _  « 

For  good  command  following,  jf($)  •  jr(s),  and  tot 
disturbance  rejection,  the  effect  of  £(s)  must  be 
kept  small.  If  the  minimum  singular  value  of 
C(s)IC(s)  Is  large  with  respect  to  unity  at  frequen- 
?las** below  crossover,  both  of  these  requirements 
can  be  met.  Likewise,  for  frequencies  above 


x^(t)  "  A  x(t)  +  £  u(t) 


(6) 


This  description  Is  modified  to  reflect  the 
(fictitious)  process  and  measurement  noises 


x(t)  -  A  x(t)  t-  L  £(t)  (7) 

X(t)  “  C  x(t)  +  £(t) , 
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what*: 


i<t)  -  pcoeaaa 
MCrlx 


with  _I  tnt«Mlty 


£(c)  ■  MaauraMnc  whlta  aolta  with  lacav 
”  alty  MCrtx.  ~ 

Th«  dtalgn  paraaatan  p  and  JL  aca  uaad  to  produea 
tha  daalrad  loop  ahapaa  of*~tha  traaafar  aatrtx 
<V>(a)  whara 

^,(a)  -  Clai  -  a]"*H  <9) 

H  -  (1//5)EC*  ,  (10) 

and  t  la  tha  aolutlon  to  tha  Plltar  Algobralc 
Xlce^tl  Aquation  (FARE) 

0  -  W  £A'  U,*  -  (l//p)W’a  .  (11) 

For  a  apaclflc  valua  of  (i,  tha  traaafar  aatrlx 
,gip(a)  can  ba  approxlaatad  qulta  raadllp.  Slnca  at 
high  and  low  fraquanclaa  a  ~  Ju,  wo  hava  tha 
approxlnatlon  (2),  (3) 

j^(a)  -  tor 

GyQ^(a)  -  C[ai  -  a]"*!,  Chen  (12) 

®l[^p*)]  *  (l/''t*)oi[^Q(^(a)]  , 

then  Che  jL  natrlx  (dealgn  paroaacer)  can  ba  choaan 
In  a  way  To  produce  the  daalrad  loop  ahapaa  and  p 
can  chan  ba  uaad  to  adjuat  the  alngular  valuaa  up 
or  down  to  Mat  the  required  croaaovar  frequancr 
apaclflcatlona. 

Aa  long  aa  la  acablllsabla  and  [a,^]  la 

daceccable,  Chen  any  choice  of  p  and  t  will  provide 
the  following  guaranteed  properclea  for  the  "target 
dealgn’ ^p(a): 

1.  cloned  loop  atable 

2.  robuatneaa  and  perfomance  guaranteen: 

"■InLL^SKpt*)]  >  1 
».ln[i+  C^p‘‘(*)]  >  Vi 

3.  Infinite  upward  gain  aargln 

4.  6  dB  downward  gain  aargln 

5.  tbO*  phaae  oarglna 

In  our  dealgn,  we  decided  to  aelect  ^  ao  chat 
all  alngular  valuaa  of  Che  target  dealgn  ^p(J<d) 
are  Identical  at  lew  and  high  frequenclea,  chua 
hoping  Chat  all  alngular  valuea  will  croaa  at  about 
Che  aoae  frequency.  To  find  Che  nuaerlcal  value  of 
L  Chat  will  aceoapllah  thla  we  proceed  aa  followa. 

Recall  that  C(a)  -  Cp(a)C,(a),  and  define 
C(a)  -  C[al_  -  a]“Ib,  where 


0  0  1 

A  -  “  ”  C  ■  [O  &]• 

•1  1 

al-A  ■  ~  ,  and 

-Bp  al-^j 

1  r  I/a  o' 

[al-A]  -  .  “  .  . 

Ac  low  fraquanclaa,  a^-^  ■ 

[•i."*p]“*  "  Since  ^  baa  dlaclncC  and  non- 

aero  elganvaluea,  exleta,  Ha  now  partition 

Che  L  natrlx  Into  ^  and  where  ^  will  be 
eele^ed  for  aacchlng  Che  alngular  valuea  at  low 
fraquency  and  ^  will  be  aelected  for  natchlng  the 
alngular  valuea  at  high  frequency, 

Fomlng  ^Q|_(*)  for  low  frequenclea, 

-  C[al-ir*l. 


w> 


[2  Cp] 


"  lf~‘ 

p/a  \ 


^Bp]a/»  •  *12 


It  la  now  aeen  that  Che  alngular  valuea  can  be 
nacched  at  low  frequenclea  If  we  aelect  Che  natrlx 
^  aa  followa: 


^  At  high  frequenclea,  a_I-^  .  ,i, 

[•i-Ap]  *  !/*•  forming  GyoL(»)  for  high  frequen¬ 

clea. 


w> 


[o  Cp] 


Bp/a‘  I/a 


■  3>Bplil/a  Cp^/a  , 

The  alngular  valuea  can  now  be  natched  at  high  fre¬ 
quenclea  If  we  aelect  ^2  **  followa: 

1:2  -  Cp'(3,Cp')"*,  (15) 

2 

alnce  aa  a  ->  >,  1/a  >  1/a  ,  and  Che  aecond  ten 
donlnatea  the  naxlnun  alngular  valuea. 

The  above  nechod  for  conetrucClng  Che  ^  natrlx 
provldea  the  dealgner  with  a  guarantee  of  Identical 
behavior  of  Che  Kalman  Filter  loop  alngular  valuea 
at  both  high  and  low  frequenclea.  However,  thla 
nechod  doee  not  provide  an  opportunity  Co  directly 
control  Che  ahape  of  Che  alngular  valuea  at 
croeaover. 

Once  Che  L  natrlx  la  decenlned,  the  dealgn 
paraneter  p  la'~uaed  to  move  the  alngular  value 
ploca  up  or  do«rn  to  obtain  the  dealred  croaaover 


m 
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fr«4u*ae7«  Than  wa  can  aolva  cha  FAKE  and  caleu- 
lata  0|fa(a).  Tha  final  aalua  of  it  Chat  wa  uaad  for 
Cha  aodal  during  Cha  Kalaan  Fllcar  daalgn  proeaaa 
waa  It  ■  A. 

Tha  fourch  and  final  acap  of  Cha  LQG/LTK  daalgn 
proeaaa  Involaaa  cha  'racoaary*  of  Cha  cargac 
daalgn  GyaCa)  by  cha  eonpanaacad  plane  cranafar 
aacrls  ^^)|C(a)«  Thla  la  dona  by  aolvlng  cha 
Concrol 'ttgArale  Kleeaccl  Equaclon  (CAKE) 

0  -  -EA  -  A'JC  -  ^'C  ♦  KBB'K,  for  q  >  0.  (16) 

Dalng  cha  daalgn  paranacar  q,  wa  aolva  for  cha  JC 
aacrlx,  and  dacamlna  cha  concrol  gain  aacrlx 

C  -  2'K  .  (17) 

For  a  valid  aoluclon  of  Cha  CARE,  chraa  condlclona 

ara  nacaaaary: 

1«  [a,^]  nuac  ba  aCabllltable, 

2.  [a,c]  Buac  ba  dacaccabla,  and 

3.  Tha  noalnal  daalgn  plant  nuac  not  have 
non-nlnlaua  phaaa  zaroa. 

Providing  cha  plane  la  alnlaua  phaaa  (3),  cha 
alngular  valuaa  of  ^(a)K(a)  convarga  Co  cha  alngu- 
lar  valuaa  of  ^pla)* aa~Cha  daalgn  paraaacar  q 
Abova  croaaovar  fraquanclaa,  additional  rolloff  la 
producad  by  cha  raeovary  phaaa,  which  furchar 
anhancaa  cha  high  fraquancy  robuacnaaa  eharae- 
CarlaClea.  Aa  a  raaule,  Cha  loop  ahapa  of  G|cp(a) 
la  approxlnaCaly  raeovarad,  and  Cha  taaulclng 
eoncrollar  will  hava  cha  daalrad  parfonaanea 
eharaecarlaclea.  A  valua  of  q  ■  1000  waa  uaad  for 
our  nodal.  Tha  nlnlaun  and  naxlnua  croaaovar  fra¬ 
quanclaa  Cumad  ouc  to  ba  0.2  rad/aac  and  O.S 
rad/aac,  raapacclvaly . 

When  calculaCad  ualng  Cha  above  procedure,  Cha 
Filter  gain  nacrlx  ]l  and  the  Concrol  gain  nacrlx  0 
define  a  apeclal  type  of  conpensacor  known  aa  an 
LQC/LTR  CoapenaaCor,  dealgnaCad  aa£(a).  Thla  coia- 
panaacor  dlffera  froa  other  LQG  crapanaacors  only 
In  cha  aanner  In  which  £  and  ^  are  calculated.  The 
acata  apace  daacrlpCloiT of  Cha  iqc/LTR  coapanaacor 
la 

z(t)  -  (A  -  BG  -  HC)  z_(t)  -  H  e(t)  (18) 

u(t)  -  -  C  z(t)  ,  (19) 

and  la  ahown  plctorlally  In  Figure  6. 


Figure  6  Scaca  Space  Deacrlpclon  of  cha 
LQC/LTR  Daalgn 


Deacrlpclon  of  Cha  Active  Roll  Control  Daalgn 

In  Che  active  roll  control  daalgn  we  uaad  all 
four  available  controla,  naaely  Che  (aealed)  bow, 
rudder,  and  dlffaranclal  acem  plane  deflaeclona. 
Theae  four  Indapendanc  concrol  varlablea  can  be 
uaad  to  connand  Che  four  (aealed)  oucpuca,  l.a., 
roll  angle,  plceh  angle,  yaw  race,  and  depch  race. 

Tha  ouccoaa  of  Che  L()G/LTR  daalgn  for  Chla 
ayacaa  la  aiuarlzed  In  Che  alngular  value  ploca  of 
Figure  7.  Flgurea  7(c)  and  7(d)  ahow  how  cha  ape- 
clflc  choice  of  Che  daalgn  nacrlx  ^  InpacCa  Che 
'target  daalgn*;  noclea  Chat  ac  both  high  and  low 
fraquanclaa  all  alngular  valuea  ara  eaaantlally  tha 
aaaa,  and  chat  tha  croaaovar  frequanclea  ara  In  the 
aane  ballpark,  aa  daalrad.  Conparlaon  of  Flgurea 
7(d)  and  7(a)  llluacracaa  how  cha  LTR  worka.  The 
singular  values  of  cha  LQG/LTR  loop  transfer  func¬ 
tion  nacrlx  are  very  near  Chose  of  the  Kalnan 
filter  loop  until  about  1  rad/sec.  After  chat  fre¬ 
quency  the  singular  valuaa  roll-off  aC  -40  dB/dac 
providing  extra  robuatnaas  to  hlgh-fraquancy  nodal 
arrora  and  sanaor  nolsaa. 

The  excellent  connand-f allowing  properClea  of 
this  design  becone  apparent  fron  Figure  7(f)  which 
shows  the  singular  values  of  the  closed  loop 
sysCen,  fron  connand  Inputs  Co  outputs.  Notice 
that  all  sinusoidal  connand  Inputs  at  frequencies 
below  0.1  rad/sec  ara  followed  with  nlnlnal 
cracking  error. 

Description  of  tha  Dealen  Without  Active  Roll 
(antrol 


If  Che  stem  planes  ara  slaved  togachar,  chan 
wa  only  have  Chraa  Independent  controls,  l.a.,  bow, 
rudder,  and  stem  concrol  surface  daflactlons. 
Thus  we  can  Independently  concrol  three  outputs, 
l.a.,  pitch  angle,  yaw  rata,  and  depch  rate.  We 
cannot  concrol  directly  Che  roll  angle,  and  wa  uusc 
accept  whatever  roll  angle  Is  generated. 

Figure  8  provides  a  suaoary  of  the  LQG/LTR 
design  for  this  3-lnpuC  3-oucput  feedback  control 
systen.  Conparlsons  of  Figures  7  and  8  shows  Chat 
Che  resulting  feedback  system  meets  essentially  Che 
same  specifications  In  Che  frequency  domain  as  Che 
4-lnpuc  4-ouCpuc  design,  except  that  we  have  lost 
our  ability  to  directly  Influence  Che  roll  angle. 
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KOHLIKEAR  SIMaATIONS 


Introduction 

Tho  eontrollor  dootgno  worn  tooted  using  both 
the  llnsor  and  non-llneor  ■ubnerslble  slaulatlons 
to  dotomlna  how  closely  the  perfomance  specifica¬ 
tions  are  net,  and  to  test  for  Instabilities  In  the 
design.  In  both  cases .  the  conpenaator  designs 
were  shown  to  be  satisfactory. 

Active  Boll  vs  Non-roll  Control  Designs 

Having  established  the  validity  of  the  conpen- 
sator  design  and  the  conpenaator  software.  It  Is 
now  necessary  to  deaonstrata  the  perfomance 
characteristics  of  the  active  roll  control  feedback 
systen  to  the  equivalently  designed  aysten  without 
roll  control  capability. 

Conparlsons  are  shown  for  two  slnulatlona.  The 
first  conparlson  la  for  a  conblned  naneuver  In 
which  step  connanda  of  1  degree  of  pitch,  O.S 
fest/second  of  depth  rate,  and  1  degree/second  of 
yaw  rata  are  provided  to  the  feedback  ayatema.  The 
second  simulation  la  for  a  connanded  yaw  rate  of  3 
degreas/aacond,  and  provides  additional  insight 
Into  the  differences  In  the  two  coopensators,  and 
Che  robustness  of  the  conpenaator  design.  The 
Interested  reader  can  find  sore  slnulatlona  In  (1) 
and  (8). 

In  both  slnulatlona  the  connanda  are  applied  as 
ranped  step  Inputs  at  C  -  S  seconds. 

Conblned  Maneuver 

This  naneuver  Is  for  step  Input  connands  of 
sero  roll  angle,  -O.S  ft/sec  In  depth  rate,  1 
deg/sec  In  yaw  rate,  and  1  degree  In  pitch. 
Referring  to  Figure  9,  we  observe  a  6Z  decrease  In 
the  forward  velocity. 

Looking  first  at  the  systen  with  roll  control. 
It  Is  observed  that  the  errors  In  roll  angle  and 
yaw  rate  are  damped  by  t  -  AO  seconds.  The  errors 
In  pitch  and  depth  rate,  however,  are  not  daoped 
until  t  -  140  seconds.  By  t  ■  200  seconds,  the 
ship  has  experienced  a  depth  rise  of  80  feet.  The 
stern  planes  are  deflected  differentially  to  coun¬ 
teract  the  roll  nonent,  with  a  steady  differential 
deflection  of  6*.  The  bow  planes  are  deflected  at 
-1.5*  to  nalncaln  the  conswnded  depth  rate,  and  the 
rudder  la  deflected  at  -2*  to  maintain  the  com¬ 
manded  yaw  race. 

Comparing  the  design  without  roll  control.  It 
Is  observed  chat  the  ship  experiences  a  snap  roll 
of  10*.  This  roll  angle  causes  a  pitch  angle  of 
-2*  which  results  In  a  large  pitch  error.  In  fact, 
at  t  ■  200  seconds,  there  Is  still  an  error  In 
pitch  of  0.5*,  or  50Z  of  the  commanded  pitch  angle. 
This  also  causes  a  -0.35  fc/sec  depth  rate  Instead 
of  Che  commanded  -.5  ft/sec.  The  net  result  of 
these  errors  Is  displayed  In  the  depth  of  the  ship. 
The  depth  rise  In  this  design  Is  25  feet.  Instead 
of  80  feet,  as  In  the  model  with  roll  control. 
Note  here,  that  a  depth  rise  of  (190  x  0.5  -)  95  ft 
Is  connanded. 


Tha  ataady  stata  stam  planaa  angla  la  -0.75*, 
which  Indicates  tha  stam  planes  are  being  used  to 
obtain  the  ordered  pitch  angle.  Because  tha  depth 
rate  Is  a  raault  of  the  conblnatlon  of  pitch  angle 
and  ahlp's  apeed,  we  observe  the  bow  planas  are 
essentially  being  used  to  attain  the  orterad  depth 
rate.  In  tha  roll  control  design,  the  ship 
obtained  the  ordered  pitch  angle  rather  quickly, 
thua,  Che  bow  planes  are  deflected  In  the  opposite 
direction  to  Unit  the  depth  rate  to  -0.5  ft/sec. 

Hard  Turning  Maneuver 

This  naneuver  la  for  a  connanded  yaw  rate  of  3 
deg/sec,  and  Is  provided  to  display  the  effects  of 
control  surface  saturation.  Referring  to  Figure 
10,  we  observe  a  drop  In  ship's  speed  of  almost 
45Z.  Looking  at  the  design  with  roll  control.  It 
Is  observed  that  the  ship  Initially  rolls  outward 
approxlnately  8*,  then  snaps  Inward  at  t  •  14 
seconds.  The  naxlnun  downward  pitch  angle  reaches 
4*  at  t  •  160  seconds,  and  starts  to  reduce  by  the 
end  of  the  run.  The  undealred  depth  loss  In  this 
case  is  184  ft.  The  stem  planes  again  deflect 
differentially  to  counteract  the  roll  nonent,  but 
now,  we  observe  the  port  stem  planes  are  deflected 
at  -3.9*  at  t  -  200  seconds  whereas  the  starboard 
stern  planes  are  deflected  at  7.8*.  This  Indicates 
Chat  Che  stern  planes,  although  deflecting  dif¬ 
ferentially  for  roll  control,  are  also  being 
deflected  to  control  Che  pitch  angle.  The  bow  pla¬ 
nes  are  deflected  at  6.25*  In  an  attenpt  Co  nlnl- 
nlze  depth  race.  To  maintain  tha  ordered  yaw  rata, 
the  rudder  Is  deflected  -27*  at  Che  end  of  the  run. 

Coaparlng  the  systen  without  roll  control,  we 
observe  chat  the  ship  snap  rolls  Inboard  19*,  and 
pitch  angle  approaches  -12*.  The  stem  planes 
deflect  to  Unit  depth  rate.  The  bow  planes, 
however,  saturate  In  this  run  at  c  ■  22  seconds. 
Up  Co  this  point,  the  ship's  depth  was  maintained 
fairly  well.  As  soon  as  the  bow  planes  saturate, 
the  depth  rate  Increases,  causing  the  ship  to  lose 
depth.  This  causes  the  stern  planes  Co  deflect  In 
Che  opposite  direction  In  an  attempt  to  minimize 
pitch  angle  and  depth  rate.  Ac  C  •  25  seconds,  the 
pitch  angle  steadies,  and  starts  to  come  off.  Ac 
t  -  108  seconds,  the  depth  rate  goes  negative,  and 
it  Is  observed  the  bow  planes  come  out  of  satura¬ 
tion.  By  t  •  200  seconds,  we  observe  that  the  roll 
angle  has  been  reduced  Co  8*,  maximum  negative 
pitch  angle  Is  7*,  depth  rate  Is  significantly 
reduced,  and  none  of  the  control  surfaces  are 
saturated.  Depth  at  Che  end  of  the  run  Is  820 
feet,  which  equates  Co  an  undealred  depth  loss  of 
320  feet,  as  compared  Co  Che  roll  control  model's 
depth  loss  of  184  feet. 

The  purpose  of  this  run  was  Co  demonstrate  how 
different  the  submersible's  trajectory  Is  when  the 
control  surfaces  saturate,  and  the  stability  and 
robustness  of  Che  design.  The  numerical  results 
demonstrate  once  more  the  beneficial  consequences 
of  active  roll  control  resulting  In  significant 
reduction  to  undesirable  depth  changes. 
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OONaUSIOHS 


Hulclvtrlabl*  control  tynton  design  using  ths 
LQG/Ln  nsthodology  hss  bssn  sucesssfully  uclllcsd 
to  design  a  subasrslbls  control  syatsn  with  roll 
control  capability. 

Ths  purpose  of  this  paper  was  to  dswonstrats 
the  adwantagss  of  roll  control  on  a  aubnsrslbls.  k 
Halted  nuabsr  of  slaulatlona  were  psrforasd.  and 
ths  psrforaancs  of  ths  aubaarslbls  with  roll 
control  la  such  laprovsd  over  ths  design  without 
roll  control.  The  control  syatsa  was  designed  for 
a  subaarslble  at  high  speed,  and  we  observed  the 
control  systen  did  fairly  well,  even  for  a  45Z 
decreasa  In  forward  velocity  of  the  ship. 
Additionally,  the  control  ayatea  was  designed  using 
Che  Inertial  reference  fraae  rather  chan  the  body 
reference  fraae  of  Che  ship.  Use  of  the  fixed 
Inertial  coordinate  systen  provided  better  control 
of  Che  subasrslble  In  aaneuverlng  situations'  chan 
for  previous  designs  which  used  the  body  reference 
fraae. 


AC  this  point  It  la  laporcant  to  stress  the 

following  observations: 

•  The  perforaance  characceriaclca  of  the 
subaerslble  with  active  roll  control  are 
enhanced  considerably  over  Che  design 
without  roll  control.'  The  slaulations 
deaonstrated  considerable  depth  laprove- 
aenc,  and  less  control  surface  deflections 
and  saturation  In  severe  aaneuvers,  as 
denonscrated  In  Figure  10. 

•  This  research  denonacratas  a  technique  to 
slaulate  perforaance  characteristics  of 
’paper'  control  aysCeas  for  trade-off 
studies  for  specified  perforaance  cri¬ 
teria. 

•  The  fact  that  only  snail  perturbations  can 
be  applied  In  validation  of  the  design  Is 
not  a  llalcaclon  of  Che  control  design 
nechodology.  It  Is,  however,  a  llalcaclon 
of  Che  linear  aodel. 

•  To  deaonscrate  the  flexibility  a  controls 
engineer  has  when  using  multivariable 
control,  Che  bow/f alrwater  planes  were 
Included  In  this  thesis.  Use  of  Che  bow/ 
falrwacer  planes  at  high  speed  In  current 
subaerslblas  may  not  be  considered  prac¬ 
tical  due  to  flow  disturbances,  and  struc¬ 
tural  llalcaclona. 
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APPENDIX  A 

ORIGINAL  MATRICES  PRIOR  TO  SCALING 
X  -[uvwpqt*9] 

A  MATRIX 
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